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Abstract 

We consider two different spin asymmetries in Drell-Yan processes generated by 
the collisions of an unpolarized proton beam on transversely polarized protons: the 
muon helicity asymmetry and the left-right asymmetry. We calculate the asymme- 
tries in the framework of the QCD improved parton model, taking into account the 
parton transverse momentum and considering the first order QCD corrections. The 
muon helicity asymmetry is sensitive to the quark transversity distribution and is 
nonvanishing even at zero order. On the contrary the left-right asymmetry vanishes 
at zero order but not at first order in the QCD coupling constant, as a result of gluon 
contribution. 



PACS numbers: 13.85.Qk, 13.88.+e 



1 Introduction 



A major problem concerning the spin structure of the proton is the possibility of 
inferring the quark transversity distribution [1-5] from data [6, 7]. As pointed out by 
various authors [8- 12], it is quite difficult to realize an experiment that is definitely 
sensitive to this quantity, which we shall denote by h{(x), f being the quark flavors. 
Indeed 

i) e — p deep inelastic scattering yields terms which are proportional to h[ but 
suppressed like rrif/Q, where mj is the current quark mass and Q the mass of the 
virtual photon. 

ii) The double transverse spin asymmetry |5], in Drell-Yan processes is propor- 
tional to the product h{(x)h 1 (x')^l\, which is probably very small, since we believe 
that, similarly to the helicity distributions, |/ix(a:)| << \h{(x)\. 

in) Other suggested |9], |14]] or planned |T5| experiments are rather complicated and 
may present serious mishaps. 

Up to now we can make only a rough evaluation of /i(|fL6l, based on the recent 
results of seminclusive n electroproduction from HERMES and SMC|| [TJ. There- 
fore doubts have been cast on the possibility of determining In the present 
situation, the best we can do is to search for all possible spin asymmetries, which 
are related to the transversity distribution, and to try to extract as much as possible 
information from the data, when it will be available. The aim of the present paper 
is to suggest a new experiment sensitive to h{. We propose to observe a muon pair 
produced by a collision between a transversely polarized proton and an unpolarized 
one, detecting, by means of standard experimental techniques [|17]], the muon helicity 
asymmetry A\ - that is the average longitudinal polarization - of at least one of the 
muons (say fi~). This kind of proposal is complementary to a typical double polariza- 
tion experiment (e. g. inclusive DIS with polarized proton target and electron beam) 
and is somewhat analogous to Collins' idea]l8||, which consists in detecting the final 
quark polarization through the azimuthal asymmetry of the fragmentation function. 

Spin asymmetries (especially with a single transverse spin) [19-33] are being cur- 
rently considered in the literature [B3], in view of experiments realized at FNALp4[ 
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and DESY|J, or planned at RHIC[37-40] and other facilities |p~3|, |35||. The aim of such 



proposals is to obtain the magnitude of twist-three terms, from which the spin asym- 
metries are expected to generate a nonzero contribution. The calculation of such terms 
involves technical difficulties and not all authors agree on the final results [|39|. Of par- 



ticular interest is the left-right asymmetry p5, 27], often called single transverse spin 
asymmetry, which can be determined from the Drell-Yan collision described above by 
detecting only the muon momenta. This asymmetry - which we call A2 - is nonzero 
only if we consider, at least, the first order correction in the QCD coupling constant 
g, for which the contribution of a soft gluon pole plays a major role. However two 



different calculations 25, 27 do not lead to the same formula. 



We calculate asymmetries A\ and A2 at tree level, assuming the QCD parton 



model [[40] (see also|22[), and successively inserting "real" gluon corrections, for which 
we adopt an axial gauge |4"1|, |4"2]| , allowing us to develop a formalism as close as pos- 
sible to the parton model. In calculating the gluon corrections, we introduce the 
quark-quark correlation functions, typical nonperturbative quantities. Our approach 
is different from the one by other authors 0, [25|, |27| , who adopt the formalism of 
quantum field theory. However we match the two ways, obtaining a precise definition 
of the correlation function in terms of destruction and creation operators. This sug- 
gests a general procedure for calculating higher twist contributions in parton model. 



Our results on A2 substantially agree with those of ref. p5|, however the gluonic 
pole contribution is unfounded. In our approach the correlation functions are inter- 
preted in terms of transverse quark and gluon polarization; moreover they can be 



evaluated quantitatively within a model by Qiu and Sterman [|T9|l , which we extend in 



quite a natural way to the case of polarized correlation functions. The muon helicity 
asymmetry is strictly related to the transverse momentum dependent transversity dis- 
tribution, whose integral over transverse momentum is h{. Furthermore we observe 
that neglecting transverse momentum implies the vanishing of the average helicity 
of the final muons, as a consequence of parity conservation. In this connection we 
like to stress the importance of the transverse momentum of the partons in various 
asymmetries . 

The Drell-Yan cross section is generally complicated [42-44]; it becomes consid- 



erably simpler p8|, |4"5| - and in particular it may be approximated by a convolution 
over transverse momentum ||4~6|, |47f| - if we limit ourselves to relatively large transverse 



momenta of the muon pairsj47|. This limitation is not a serious obstacle for our aims 
and it is even particularly appropriate for determining h{, as we shall see. 

In sect. 2 we define the two asymmetries we want to calculate and we give some 
general formulae for the asymmetries, the Drell-Yan cross section and the expression 
of the leptonic tensor. Sect. 3 is dedicated to the calculation of the hadronic tensor 
in the QCD parton model. In sect. 4 we write the QCD first order corrections by 
adopting an axial gauge. Moreover we impose gauge invariance and establish some 
symmetry properties for the correlation functions. Then we illustrate the properties of 
the perturbatively calculable ("hard") coefficients; lastly we perform the calculations. 
In sect. 5 we give the expressions of the asymmetries and evaluate their orders of 
magnitude. In sect. 6 we draw some short conclusions. 



2 General Formulae 

In Drell-Yan processes generated by the collisions between an unpolarized proton 
beam and a transversely polarized one, p^p — > fi + fi~X, two different kinds of asym- 
metries can be defined, i.e., the muon helicity asymmetry, if the polarization of one 
final muon is detected, and the left-right asymmetry, if only final momenta of the 
muons are determined. 

i) The muon helicity asymmetry is defined as 

da+ — da- , N 



where da± is the inclusive Drell-Yan differential cross section with the final fi~ having 
positive (negative) helicity. 

ii) The left-right asymmetry is 

^ = d(J r ~ doi ^ 
da r + doi ' 

where da r ^ is the differential cross section with a final /i~ at the right (left) of the 
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plane determined by the momentum and by the spin of the polarized proton. Indeed, 
if the muon helicities are not detected, the only way of constructing an asymmetric 
term (i.e., containing the e tensor) is 



A 2 oc S ■ p, (3) 
where p is the pT four-momentum and 

S a = e a {3 lS n^rP,S 5 , (4) 

S being the Pauli-Lubanski four-vector of the polarized proton. Here we also have 
set 

and Pi and P 2 are the four-momenta of the two protons, whose spatial parts are Pi 
and P2. Proton 1 is polarized, whereas 2 is not. In the laboratory frame, where Pi = 
-P2, eq. (|3|) implies that the number of p~ falling at the right of the plane determined 
by the spin and momentum of proton 1 is different from those which occur at the left 
of that plane. 

The Drell-Yan cross section reads 



da = T7r7= —L flv H^ dT, (6) 



1 e 4 

where dT is the phase space element, whose expression is 



dT = WT2^P S (P 2 ) *(Po) d 'P $(f) Ofo) S\P +P-Q)- (7) 
(27r) i 

Here p and p are the four-momenta of p + and yT respectively, q the four-momentum 
of the virtual photon, Q its effective mass, P = \P±\ = \P 2 \ the modulus of the 
momentum of each proton in the laboratory frame, s the overall energy square, L^ u 
the leptonic tensor and H^ u the hadronic tensor. 

The leptonic tensor is, in the massless approximation, 
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V = \Tr\j(l ±75)7^ 7,], (8) 
± being the sign of the the /i~ helicity. Then we can write 

L^ u = S 1 ^ ± iA^, (9) 

where 5' and i-Aj^, are, respectively, the symmetric (real) part and the antisymmetric 
(imaginary) part of the leptonic tensor: 

S U = Pi&v+PuPv- 9»uP-p, (10) 

A\ v = e^uV u f- (11) 
Concerning the hadronic tensor, we set 

H, u = Hj8 + H'^ (12) 

where Hj®J and H' are, respectively, the QCD zero order and first order approxima- 
tions, which will be calculated in the next two sections. 

3 Drell-Yan Hadronic Tensor in the QCD 
Zero Order Parton Model 

Now we write the hadronic tensor for the high energy Drell-Yan process. According 
to the QCD parton model 1 4(J, we take into account the transverse momentum of the 



quark inside the hadrons. Moreover, as already explained in the introduction, we 
limit ourselves to time-like photons with transverse momenta |q_i_| of order 1 GeV 
and with Q » |q_i_|, where Q 2 = q 2 and q is the four- momentum of the pair. Then 
the generalized factorization theorem ||19||, in the covariant formalism [48|, yields 



Hf) = E (13) 
where / are the three light flavors (u, d, s), e\ = 2/3, e 2 = e 3 = -1/3, and 
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K = JdT q Y: [qk(Pi)VT?(P2Kl 2 (puP2; S) + (l~ 2)]. (14) 



Ti,T 2 



Here q T (g4) (/ = 1,2) are the probability density functions of finding a quark (an- 
tiquark) in a pure spin state whose third component along the proton spin is TJ. 
Moreover dT q is the analog of the element of phase space of muons (see eq. (0)): 



dr «=7^3I dW(P?)0(Pio) d 4 p 2 5(^)%2o) 5 4 ( Pl +p 2 -g) (15) 
and are the four-momenta of the two active partons. These are taken on shell and 
massless, which is a good approximation for the values of |q_i_| considered. Lastly 

h%? = ^Tr(p T ^p T * lu ), (16) 



where the factor 1/3 comes from color averaging in q-q annihilation, T\2 = (Ti,T; 



2) 



and p is the spin density matrix P9 



P Tl (p Tl ) = \m + 2T ll5 (± m + f) ± )]. (17) 

2T^?7|| and 2Tir)± are respectively the helicity and the transverse Pauli-Lubanski four- 
vector of the active partons. We have r)\\ = Sn^, where n/ = p«/|p/| and S and pi are, 
respectively, the space components of S and of pi in the laboratory frame. Moreover 

= (0, S — 7711 iii) (I — 1; 2) and 7711, a Lorentz scalar in the limit of zero quark mass, 
can be defined covariantly |)(J . 

Carrying on the integration (|14D over the time and longitudinal components of p±, 
and adopting the light cone formalism, we get, in the limit of high Q 2 and |p;j_| << 
Q(l = l, 2), 



E U = T^7V2 / dVi E ^(^i, Pu)^(^2, p 2 i)^(^i, ^ 2 ; 5) + (1 ^ 2)]. (18) 

47T J T1T2 

Here = (go ± q\\)/\fs are the longitudinal fractional momenta of the two active 
partons, go and q\\ are respectively the time and longitudinal component of q, and 

P2± = q± - pi±. (19) 



Since the proton 2 is unpolarized, we set 



1(T 2 =l/2) - S(T a =-l/2) > Q(T 2 =l/2) ~ 9(T 2 =-l/2)- ( 20 ) 

Inserting eqs. ( |IED to ( j2"0| ) into eq. fll^D, and taking into account the relation 

Q 2 = A Xl x 2 P\ (21) 

we get 

= + (22) 

where S* and are respectively the symmetric and antisymmetric part of the 
hadronic tensor: 



*~ - W a - SQf - (24) 



where 



Q f = J(f Pl± [q[(x 1 ,pl ± )q f 2 (x 2 ,pl ± ) + (1^2)}, (25) 
6Q f = J d 2 p 1± —^^[Sq{(xi,p 1 ±)q((x 2 ,pl ± ) - Sq f 1 (x 1 ,p 1A _)q^(x 2 ,pl ± )). (26) 



Here we have introduced the distribution functions: 



1/2 1/2 

4= E s<d= E ar^, (27) 



T,=-l/2 T t =-l/2 



and similarly for the antiquarks. Furthermore 



■V = ni M n 2i/ + r^r^ - g^ u , = e afl p u n"n%. (28) 

Such tensors fulfil gauge invariance up to twist-4 terms. The distributions 5qf are 
related to the transversity functions h{: 



h{(x) = I d 2 p ± 5q f (x,p ± ), (29) 

a similar relation holding true for antiquarks. 

Now we discuss a symmetry property of the spin density functions 8q* . Invariance 
of strong interactions under parity, time reversal and rotations (in particular rotations 
of 77 around the proton momentum) imply 



Sq f (x, p±) = Sq f (x, -p±). (30) 

However, if we take into account the initial state interactions and soft gluon ex- 
change between the protons, we have to introduce also the effective, T-odd, density 
functions]^], which do not fulfil eq. (|3~0"D, since time reversal invariance does not 



apply trivially to them. This fact has no particular consequences on the integral (|26|) . 
which in general does not vanish for |q_i_| of order 1 GeV. However, if we integrate 
the cross section over the transverse momentum of the muon pair, the the antisym- 
metric hadronic tensor derives its contribution (if any) from the sole effective density 
functions. 

Inserting eq. ( |2"2"| ) into eq. flT3|), we get 



H$ — S$ + iA^ u , (31) 

where 

3 3 

/=i /=i 
In order to find the asymmetries A\ and A 2 (see eqs. (|l|) and @), we combine the 

hadronic tensor ([!!]) with the leptonic tensor (|9|), according to formula (|^) for the 



differential cross section. Therefore, in the zero order approximation of the QCD 
parton model, we find the following results: 

i) If the helicity of the final negative muon is detected, the leptonic tensor has 
a nonvanishing antisymmetric part, which, combined with the antisymmetric part of 
the hadronic tensor (0), gives a nonzero contribution to the muon helicity asymmetry 
A 1 . 
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if) On the contrary, if no helicity is detected, only the symmetric part of the 
leptonic tensor survives and combination with the hadronic tensor drops out A^J. 
Since (see eq. (|23|)) is spin independent, the left-right asymmetry A 2 vanishes. 
This result may be also deduced from parity and time reversal invariance. 

Concerning the muon helicity asymmetry, we are faced with two important ques- 
tions: 

1) The asymmetry A\ and the transversity function, h{, depend on the distribution 
function 5q? through integral relations, respectively eq. ( p6|) and eq. (f29|). Therefore, 
we have to prove that the function we extract from eq. contributes to h{. 

2) One may wonder where a does chiral-odd function like h{ come from in the 
process considered. Such a distribution function arises, e. g., in DY with two polarized 
proton beams, or from a mass term, which causes helicity flip. 

To answer the first question, we write the relation 

8qf = cosa (Sq R — 8q{ j ) + sina 8q H , (33) 

where 

5q f R(L) = \<q R[L) \Pl >\ 2 , (34) 

5q f R = 2Re{<pl\q R ><q L \pl >], (35) 

P T ' S foei\ 
cosa = — -p- . {6b) 

Here > and \qn(L) > are shorthand notations for a canonical state of the 
proton and for quark helicity states. Relation ([33]) follows from decomposing into 
helicity states a quark state with a given transverse momentum and spin parallel or 
antiparallel to the proton spin. The first term of that relation, which is twist 3 and 
chiral even, may be expressed as a linear combination of the usual helicity distribution 
functions. On the contrary, the second term is chiral odd and prevalently twist 2, 
since sina ~ 1. Thanks to arguments similar to those which have led to relation 
(130|) , we conclude that the first term of eq. ([33]) does not contribute to h{, as can 
be seen from eq. (p9|). The transversity function depends solely on the chiral-odd 
function 5q{j, which contributes to the muon helicity asymmetry, as can be checked 
by substituting eq. ([33]) into eq. (p6|). 
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As regards the origin of the chiral-odd function in our process, we show in appendix 



A that the antisymmetric tensor (|24j ) may be written as 



Ku = m f^W I d W E [<Zt>i> Pu)^(z 2i pl±)Sj ia pf - (1 <- 2)] . (37) 

Ti,T 2 

Here m./ and S'J 1 are respectively the mass and the Pauli-Lubanski four- vector of the 
active quark. The quark mass causes helicity flip and therefore a contribution to the 
chiral-odd distribution function. The smallness of my is compensated by Sj 1 , which 
results to be (see Appendix A) 

5? ~ 2T iPi _l • S^i. (38) 
1 rrifXi 

The effect we have just illustrated is washed out in DIS, since we have to convolute 
the elementary quark-photon cross section over only one distribution function. In 
this case also the initial state interactions, and therefore the effective distribution 
functions, are suppressed. 



4 QCD first-order contributions to spin 
asymmetries 

Now we consider the q — q annihilation amplitudes with one "real" gluon emitted 
(absorbed) by one of the colliding hadrons and absorbed (emitted) by the active 
parton of the other hadron. These interfere with the amplitudes just considered in 
the preceding section. The interference terms are twist-three and first order in g. We 
adopt the axial gauge A a ■ n 2 = 0, where A a is the gluonic field, a = 1...8 and n 2 is 
given by eq. (|). This gauge is not covariant, but allows to keep the parton model 
description; in particular, in our approximation, this gauge avoids the complication of 
"link" operators |52|, which we should introduce in a generic gauge, in order to ensure 



gauge invariance of the nonperturbative ("soft") functions involved in scattering. 

In this section we write the twist-three contributions according to the parton 
model[f46]; secondly we impose gauge invariance and deduce symmetry properties of 



the "soft" functions that we are going to introduce, the so-called correlation functions; 
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thirdly we discuss the behavior of the fermionic propagator, which appears in the 
"hard" coefficient; lastly we perform calculations. 



4.1 Parton Model Approach 

The twist-three contribution to the hadronic tensor reads, in tree approximation, 

H^ = gj:<? f {H'Jf-H't*), (39) 

where 

H'^^H'^, H^ = j:H'^. (40) 

i=i i=i 

Each term in eqs. Q4TD corresponds to a different graph: 

i) H'J 1 ) a ' 1 (/ = 1, 2) refer to a gluon emitted (or absorbed) by proton 1 - which is 
polarized - and absorbed (or emitted) by the active antiquark (I = 1) or quark (I = 
2) of proton 2. 

ii) H'J^' 1 refer instead to a gluon emitted (or absorbed) by proton 2 and absorbed 
(or emitted) by the active parton of proton 1. Analytically we have 

H l f f ,1= d ^i 4 1 V^ Xl ' Pl±;x/ i' p/ ^)^ 2 2 ( x 2,P2±)^ 12 (a;i,a; / i,a;2), (41) 

T lt T{,T 2 

H ' f !f' 2 =[ d ^i Yl ^^(^IjPll;^ P / i±)?4 2 ( a; 2,P2±)^ 2 {xi,x' 1 ,x 2 ), (42) 

while H'^' 1 are obtained from ( f4T|) and ( |42| ) by substituting (1 «-> 2) and // «-> v. 
c^ T , ic^rp,) (/ = 1, 2) is the correlation function between a quark (antiquark) of 
four-momentum p\ and spin T\ and another of four-momentum p\ and spin T{. dQi is 
the phase-space element, i. e., 



- m^ dWd%/ i' (43) 

where dT q is given by formula (|15l), k\ is the four-momentum of the absorbed (or 
emitted) gluon and Z[ = x\ — xi. The arrow denotes integration over p i0 , pm, k i0 , km. 
Lastly in the axial gauge A a ■ n 2 = we have 
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C 12 = l^J—Tr [k Ti ' t ^ (f 2 ^p T2 + p T2 $ 2 ) lu] , (44) 



3p 2 2 + te 
where 

Ti2 = (Ti,T(,T 2 ), p 2=P2 -k u (45) 

e M is the polarization four-vector of the gluon and k t ' t ' the spin correlation matrix 
between a quark of spin T and another with spin T'. /i„„ is obtained from eq. (fHj) 
by changing the quark matrices with the antiquark matrices and vice- versa, i. e., 
k Ti ' t i — > k Ti,t i and p T2 — > p T2 . As shown in Appendix B, 

« T > T '=v(t)<' T ', if^=m^ t (46) 



where 



3j 3/ , / \ 1 I ^ / , \ 

VKO = -*=*p (47) 



<' T = ^ = ^(1 + 2^), <'" T «'" T ) = ^75(^T2.T). (48) 



Since we are considering the twist-three contribution, we neglect the transverse mo- 
mentum in the numerators of h^ v and h . In this approximation the spin density 
matrix reads 



P T = P T = ^(i + 2T 7 5^). (49) 

In the next subsection we shall show that the sum H^J + H' - where Hf£j) is given 
by eq. (p~3|) and H' by eq. (|39|)) - is gauge invariant, although the single terms are 
not. 

Since the proton 1 is transversely polarized, quark-gluon interactions involve lin- 
early polarized gluons^9|, either in the direction of the spin of the proton, or in the 
direction orthogonal to the spin and to the momentum of the proton. In Appendix 
B we show that 



Ti = T x : e = 2T 2 S; 



T' 2 = T 2 : e = 2T x S] (50) 



T( ^ Tj : e = -S5 Ti ,t 2 - 



(51) 



13 



Substituting eqs. (fi6]) and (pESP to (|5T| ) into eq. (f0|), and developing calculations, 

we get 



hi 12 = h'gl5 Ti , Tl + h'gl5 Tli _ Tl , (52) 

where Ti 2 = (Ti,T 2 ) and 



hf,l = \^(h)'^-^2(T 1 + T 2 ){p lfl S v +p lv S,), (53) 



tiftl = ^(h)"^-^ 2T 1 e a ^ v S a pP + 2T 2 ( Plfl S„+p lv S IM )\6 TltT2 . (54) 



Here the suffixes AT and F denote repectively spin nonflip and flip of the active 
quark. The tensors ^a^fw are ec L ua l to expressions (]53f) - fl5"4|) , except for a change 
of sign in front of the first term of eq. ([54j) . h^p^ v and h N(F)p,vi involved in the 
expressions of H'J^' 1 , are obtained from k^fo and h! ^ F ^ v by substituting (1 <->• 2) 
and //<->• za Now we establish some important relations, which considerably simplify 
our calculations. 



4.2 Correlation Functions: Normalization and 
Symmetry Properties 

In Appendix C we show that gauge invariance implies the following normalization for 
the quark correlation functions: 



^Ux^^x'^Wpi ~ h-v'i) = STr^dp.llk.llp',!) 1 ^ 

x E E < ^l4 i (p / >^( k 04 iJ (p/)l^ > • (ss) 

ij=\ c =l 

Here z, j and c are the color indices of the quarks and of the gluon, combined in 
such a way that the operator product constitutes a color singlet. c^ Ti and are the 
destruction and creation operators for the quarks, a^(k^) (m = 1,2) are destruction 
operators for gluons whose polarization four- vector is S (m = 1, corresponding to 

14 



T{ 7^ Ti), or 5 (m = 2, corresponding to T( = 7]). a c m {ki) has to be substituted by 
a^(k;) for negative values of the gluon energy. 

Invariance under parity inversion, time reversal, and rotation by it around the 
proton momentum imply, through eq. (|55|), 

c^ jTi (x[, pjj_; x h pi±) = -c^ T ,(x h -pi±; x' h -p' l± ). (56) 

To show that, we observe that the combined action of the three transformations 
leaves spin and longitudinal momenta unchanged, while inverting the arguments of 
the function and transverse momenta. Moreover time reversal produces the same 
change of phase (nn, with n integer) as the rotation. Therefore the minus sign in 
front the r.h.s. of (^) is due to the product of the intrinsic parities of the two fermions 
and of the gluon. 

Furthermore the hermiticity condition yields 

c T' T X x 'n Pz±; x h Via.) = c S J* T ,{x h p /± ; x' t , p' l± ). (57) 

Therefore, if we integrate over transverse momenta, the correlation functions are 
imaginary and antisymmetric under the simultaneous exchange of the momenta and 
spin of the two quarks. Relations similar to eqs. (|53| ) - ( |5"TD hold true for the antiquark 
correlation functions. Lastly, since the proton 2 is unpolarized, we have 



J,2 _ f,2 _/,2 _ _f,2 

"t 2 ,T^ — C -T 2 -T^ C T 2 ,T^ ~ C -T 2 -T/ ! - 



These equations, together with eqs. (^(]), imply that the terms proportional to T 2 
in eqs. (p3|) and (^4|) (or in those which are obtained by substituting (1 «-> 2)) are 
dropped after summing over the spin indices. 

The symmetry property (^6|) does not hold if we take into account the initial state 
interaction and soft gluon exchange between the protons. Indeed, analogously to the 
effective density functions, it makes sense to introduce effective, T-odd, correlation 



functions, for which time reversal invariance cannot be expressed in a trivial way ||22|| . 
Therefore such functions - whose importance will be shown in the following subsection 
- have in general a real part even after integration over transverse momentum. 
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4.3 The fermionic propagator 



In twist-3 approximation the tensor (44) - the "hard" coefficient of the hadronic 
tensor - is independent of the transverse momentum of the active parton. Therefore 
our preceding considerations imply that, aside from the effective correlation functions, 
the left-right asymmetry of the cross section integrated over the transverse momentum 
of the muon pair receives a contribution from first order corrections, only if the tensor 
fl4*4] ) has an imaginary part. This is in agreement with the observation by Boer et 
al.|25|], although we arrive at a different conclusion. We examine the question in 
detail. 

The tensor (|4]) includes the factor (2ki ■ Vi)l^p2 + ^ e )- For a massless quark 
P2 = — 2&4 • £>2j which annihilates the effects of the imaginary part of the tensor, 
unless the correlation function c^ l T , has a simple pole - a gluonic po/e|25|] - just at 
p\ = 0. If we would neglect, like Boer et al. ||25|| , the transverse momentum effects in 
the fraction above, p\ would vanish at z\ = and it would make sense to assume a 
gluonic pole just at zero momentum (a soft gluon contribution), so that the fermion 
propagator would yield an imaginary part. But, at least in the approximation of on- 
shell and massless quarks - adopted also by Boer et al. -, the transverse momentum 
has dramatic consequences on the pole, which results to be located at 



^i = X2jPP2± p k 2 1±±VA , A = (^P 2 + pL)[( P2 x-k 1± ) 2 -pLk?J. (59) 

This value is generally complex, and by no means can be approximated by Z\ ~ 0. 
Consequently a soft gluonic pole would not imply an imaginary part for the tensor 
(|44l) . On the other hand it appears quite arbitrary to assume that c^, has a pole 
located just in correspondence of the value given by the first eq. fl59|). Analogous 
considerations can be done for the pole (p 2 ) -1 , with p 1 = p\ — k 2 . Therefore the 
"hard" coefficient of the hadronic tensor has not an imaginary part. We just set the 
fractions in front of the tensors fl53|) and (|54]) equal to -1. As a consequence the twist- 
3 contribution (if any) to the left-right asymmetry of the integrated cross section is 
entirely due to the effective correlation functions. Our considerations are somewhat 



analogous to those of Anselmino et al.[22|, concerning pion inclusive production. 
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4.4 Calculation of the gluon correction 



Now we perform the calculation of the hadronic tensor fl3~9"|). Substituting eqs. ([4~l|) 
to ([43]) into eqs. fl40|), and taking into account the results (|53"D and ( p4[) and the 
considerations above, we get 



H 'f,a 



H 'f,b 



6(2vr) 5 Q 2 
V2x 2 P 
6(2vr) 5 g 2 



(s 2f ,„5 2 Q f + ia 2f ,JiQ f ). 



(60) 
(61) 



Here we have set 

~5xQ f = J d 2 p 1± [6C£ tl (x u Pi±)q( (x 2 , pl ± ) +5C f SA (x 1 ,pi±)q£(x 2 ,pl ± )], (62) 
5 2 Q f 
SiQ f 



d 2 pi±[C f S2 (x 2 ,p 2± )5q{(x 1 ,p 1± ) + C S)2 (x 2 ,p 2 ±)8q{(xi,p 1± )}, (63) 
J d 2 p 1± [SC f A1 (x 1 ,p 1± )q f 2 (x 2 ,p 2 2± ) - 5C f A1 (x 1 ,p 1± )q((x 2 ,pl ± )), (64) 



d'fu.P'^pMil.Vix) -^.hW^RiJI (65) 



and 



Moreover 

5 P2_L, 

where 



^2TC^(x 1 ,p 1± ), 

T 

J2 C S0>( X 2,P2X), 
T 

X;2rC^(a; 1 ,p LL ), 

T 

X! Ca 2 A X 2, P2±), 



(1 = 1,2). 



(66) 

(67) 
(68) 
(69) 
(70) 



C^((x h pi ± ) = RelJd 2 k l± J 



dzi 

Zl 

dzi 



V>(*i) , (71) 



C$' r ((a:j, pi ± ) = Jm j y d 2 A; u y — [c£'_ T (x z , p /± ; x'„ p{ ± )] ^(t/) | (72) 
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and I = 1,2. Analogous expressions hold for the barred quantities. Here z and 
are defined in the expression fl43|) of the phase-space element, while t is given by 
the first eq. (|47|) . The QCD first order correction H' is obtained substituting eqs. 



(|60| ) and ( |6T1) into eq. (|39D ). Three remarks are in order. 



The functions Cg, 5Cg, and SC A have the dimensions of a mass, as can be 



checked from eq. fl55[). 

Integrating the cross section over the transverse momentum of the muon pair 
amounts to independently integrating over the quark transverse momentum the 
"soft" functions involved in the hadronic tensor. In particular, as stressed in the 
previous subsection, the integrals of the functions 5iQf (/ = 1, 2, see eqs. (|62j) 
and (^)) derive their contributions from the sole effective correlation functions. 

Taking into account eqs. (pE?D, the integrals ( |7T|) and ( [72] ) can be split into two 
terms, corresponding to 

-[cf(x, x + z) + c~f(x, x — z)\ and — [c^(x,x + z) — c^(x,x — z)], (73) 

where = c^ T /(l — t 2 ) -1 / 2 , = t/zc? and t is given by the first eq. (|47|). If 
we assume c/(x, x') = Cq(x,x')5(x — x'), where Cq is a smooth function of its 
arguments, the second term ([73]) corresponds to the derivative term by Qiu and 
Stermanpl, M. (See ref.[Bl for details). 



5 Asymmetries 

In QCD first order approximation the hadronic tensor H^ v consists, according to 
eq. ([12]), of the sum of two terms, which have been calculated in the two previous 
sections (formulae (|3lD and (p9|) ). Substituting the expressions of H^ u and of the 
leptonic tensor ([|) into the differential cross section (||), and taking into account the 
formulae (|1|) and (0) of the two asymmetries, we get 

A 1 = 4° ) +4 1) , (74) 



18 



where 



and 



(o) 8^x^E f e}5Qf cos9 

1 ~ Q E f e 2 f Qfl + cos*e { } 



Ai) _ ~9 £/ e)(xi5iQ f + x 2 5 2 Q f ) 2sin9sin<j> 

1 = (2vr) 3 v /^Q £/ e}Q/ 1 + cosW • ( j 

Moreover 

_ g J2j ej(xi5iQ f + x 2 5 2 Q f ) sin29cos(p 

2 = (2n) 3 y /x^Q Z f e}Qf 1 + cos 2 6 ' ( } 
Here 9 is the polar angle and 4> the azimuthal angle of the negative muon: we 

have assumed a reference frame in the center-of-mass system of the muon pair, whose 
z-axis is taken along the momentum of the polarized proton, while the x-axis is along 
the space component of S. Furthermore Qf ,8Qf ,6iQf ,5iQf (I = 1,2) are given, 
respectively, by eqs. fl25|) , fl26|) and fl62| ) to (|65|). 



The expression of A 2 turns out to coincide with the one by Boer et al.[^5| [52 
provided 

- we integrate the cross section over the transverse momentum of the virtual 
photon; 

- we take into account the following notation differences: q(x) — > fi(x), (f) — > 
$c — — 



we identify 



M P h f (x) = -9j^q J d 2 P± C f s (x, P± ), (78) 
MpfUx) = -97^q J d 2 p ± 5C f s (x, p.0, (79) 



(2vr) 3 g 

where Mp is the proton rest mass. 



5.1 Zero order approximation 

As we have seen, the left-right asymmetry at zero order vanishes. In the same approx- 
imation the muon helicity asymmetry is sensitive to the p_i_-dependent transversity 
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distributions 5qf(xi,pi±) and 6qf(xx, puj, which can be used for calculating the 
transversity distributions h[ and h x (see eq. ( p9|) ), with / = 1,2,3. Eq. (|75|) is a 
linear combination of the six unknown functions. In order to extract them, we need 
other, independent combinations. For example, considering Drell-Yan events from 
collisions between a polarized proton beam and, say, a pion, we get an expression of 
the asymmetry analogous to Ai (eq. (0)), where qf and qf are replaced respectively 
by the quark and antiquark density functions inside the pion. It is worth observing 
that Drell-Yan offers, at least in principle, a variety of independent combinations. For 
example, we could consider collisions of transversely polarized protons on unpolarized 
protons, antiprotons, positive and negative pions and kaons. In this way we could 
obtain six independent asymmetries similar to formula (|75|), from which (or from 
part of them) we could extract the unknown distribution functions by a fit, taking 
into account symmetry properties like isospin invariance, or general constraints like 
Soffer's inequality [^3| . 

In this connection it is worth recalling that in ref . ]54]] it was proposed, quite anal- 
ogously to us, to extract the quark helicity distributions from Drell-Yan produced in 
scattering on longitudinally polarized proton target of beams of pions and unpolar- 
ized protons. They obtain an asymmetry formula similar to Ai (eq. (|75|)), although 
this is a twist-two and not a twist-three effect. 

It is known that our procedure of convoluting the elementary cross section over 
the transverse momentum of the active quarks is a good approximation only for 
sufficiently large transverse momenta^]. But this is not a severe constraint, since, 
as we have seen, transverse momentum is essential for exhibiting helicity asymmetry. 
Eq. (|75|) suggests that directions not too far from the forward and backward one are 
the most proper. 

In order to estimate the order of magnitude of the asymmetry Ai, we take into 
account the HERMES results!!, i. e., \5q f \/q f ~ \5q f \/q f ~ (50 ±30)%. Moreover 



eq. ( ff5|) implies A\ to vanish for x\ = x 2 - Therefore, in order to maximize the 
asymmetry, we should take x\ as different as possible from X2, without making 
and Qf too small. For Q 2 of order 10 GeV 2 and s ~ 100 GeV 2 , a good choice, 
consistent with eq. (|2l|), is X\ ~ 0.5, x 2 ~ 0.2, or, viceversa, x 2 ~ 0.5, X\ ~ 0.2. 
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The evaluation of the asymmetry is particularly complicated, because the integrand 



at the r.h.s. of eq. fl26|) is partly positive and partly negative. We just give an upper 
limit: since |pj_| ~ 0.5 GeV, under the conditions illustrated above we estimate \Ai\ 
< (10 ±6)%. 

5.2 First order corrections 

In order to find the order of magnitude of the first order correction to the two asym- 



metries, we generalize the Qiu-StermanfT9| guess: 



c s,t( x > P±) = k sQt( x , P±), c a,t( x , P-l) = K aQt( x , P±), (8°) 

where Cg T and C AT are given by eqs. ( |7T| ) and ([72] ) and Kg and K A are constants 
and independent of T. Then eqs. flB7|) to ( f70| ) imply 

Cfe p ± ) = tf s /(ar, pi), <fCfe p ± ) = K s 8q f (x, p ± ), (81) 

analogous relations holding true for C A and 8C A . Two remarks are in order. 

• The generalization (|80|) of the Qiu-Sterman model is quite natural and imme- 
diate in our approach, whereas it would be rather complicated in the formalism 
of quantum field theory. 

• To be precise, we have generalized one of the two guesses proposed by Qiu and 



Sterman, the other one being obtained by multiplying the r. h. s. of eqs. (pi 



by x. The two models could be distinguished by examining data at very low x. 



According to refs.[19, 27, M], we may set g\K s \ = 0.08(27r) 3 v / 2M P . On the other 
hand, since the sums ( p7[ ) and QSg|) consist of four terms, whereas the sums ( |B§D and 
([70|) consist of two terms, in order to evaluate the order of magnitude, it seems not 
completely unreasonable to guess that \K A \ ~ l/2\Ks\. 

We fix some parameters as before: Q 2 ~ 10 GeV 2 , s ~ 100 GeV 2 , x\ = 0.5 and 
X2 = 0.2 (or viceversa). Moreover we take into account once more the above men- 
tioned results by HERMES. Lastly, as regards the angular dependence, eqs. (Jf/j) and 
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(|76| ) suggest that the most favourable conditions for detecting the left-right asymme- 
try and the first order correction to the muon helicity asymmetry are, respectively, 
6 ~ ~ tt/2 and 9 ~ tt/4, 3vr/4, = 0, vr. We find \A 2 \ ~ (2 ± 1)% and (A^l ~ 
(1±0.6)%. The uncertainty on these two asymmetries could be even larger, if we take 
into account that the coupling constant g is poorly known |2T|j. It is worth noticing 
that the first order correction to the muon helicity asymmetry is much smaller than 
, moreover it can be disentangled, since it exhibits a completely different angular 
dependence. 



6 Conclusions 

We have considered the muon helicity asymmetry (Ai) and the left-right asymmetry 
(A 2 ) in the Drell-Yan process pp^ —>■ fi + fi~X. The two quantities are measurable, 



although for A\ some special care is required Wt\. Let us recall the main results. 



1) We have calculated the two asymmetries in tree approximation by means of 



the QCD improved parton model [46], P0[ , which amounts to replacing transverse mo- 
mentum by a combination of momentum and gluon field, as dictated by the covariant 
derivative D a (see formula ( |C.3| )). Indeed both asymmetries vanish at the leading 
twist and turn out to be generated, at twist-three, by the transverse components of 
the covariant derivative, i. e., 



D a S a -> A u D a S a -> A 2 . 

In particular A± receives its main contribution from the first term of the covariant 
derivative, that is, from parton transverse momentum, while A 2 gets a contribution 
from the sole second term, corresponding to one gluon corrections. 

2) We have performed the calculation in the framework of parton model, assuming 
an axial gauge and introducing the correlation functions between two partons of dif- 
ferent momenta in the hadron. We have matched this simple and intuitive formalism 
to the quantum field theory approach, obtaining, as a result, condition fl55|), which 
uniquely fixes the normalization of the correlation functions and simultaneously guar- 
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antees gauge invariance. This procedure is of quite general validity for twist-three 
terms, owing to the local character of the interaction. 

3) If we integrate the cross section over the transverse momentum of the muon 
pair, only the effective density and correlation functions contribute to A\ and A%. 
Therefore it appears useful to consider these kinds of asymmetries as well as those 
from the differential cross section. 

4) Our approach leads to a simple physical interpretation of the "soft" functions 
C* and 5C? involved in the left-right asymmetry, and to a natural generalization of 



the Qiu-Sterman[|r9| assumption, allowing for a quantitative evaluation of A 2 . As 
a result, for standard values of Q 2 and s, this asymmetry is estimated to be a few 
percent. 

5) Under the same conditions the muon helicity asymmetry is probably a bit 
larger. This asymmetry is sensitive to the transverse momentum dependent transver- 
sity distributions, 5qf and 5qf , which are related to the usual transversity distribu- 
tions through eq. (|29|). One possible advantage of our method is that one can obtain 
independent combinations of such functions by performing various scattering experi- 
ments with beams of unpolarized protons, pions and kaons on a transversely polarized 
proton target. 

Acknowledgments: The author is deeply indebted to his friend prof. J. Soffer 
for useful and stimulating discussions on this work. Furthermore he is grateful to his 
friend prof. C. M. Becchi for important suggestions. 



Appendix A 

Here we calculate the antisymmetric part of the hadronic tensor, by inserting the 
mass of the quark in the density matrix. We show that it reduces to the expression 
(24) in the limit of zero mass, the smallness of the mass being compensated by the 
"good" component of the Pauli-Lubanski (PL) four- vector of the active quark. The 



23 



hadronic tensor reads 



E U = -^02 1 E [4 1 (^i ) P 1 ±)94(^,p 2 ±)4 2 (xi,x 2 ; S) + (1 «-> 2)], (A.l) 



4vr 2 Q 

where 
and 



Ti,T 2 



/# a (z ls ^2; 5) = iTr(p Tl 7 ,p T2 7,) (A.2) 



p T (p T )= 1 -^±m f )(l + 1 ^J). (A.3) 
Here p = p\ or p 2 and the + and - sign refer, respectively, to the quark and to the 



antiquark; moreover Sj is the PL four-vector of the quark. We have 



S% = 2TL%. (A.4) 

Here £ = (0, S) is the PL four-vector of the proton, in the frame where it is at rest. 
L is the matrix of the boost which takes the quark from rest to a four-momentum 



p = (yrriq + p 2 , p), where p = (p_i_, xP) and P is the proton momentum. A standard 
calculation shows that 

IT 

P±-S(l, 0,0,1). (A.5) 

1 rrifX 



On the other hand the antisymmetric part of the tensor ( |A.1|) turns out to be 

A U = J d2 Pi± £ alulvhix^P^hi^P^) ~ i 1 2 )]> ( A - 6 ) 



4vr 2 g 

where 



Ti,T 2 



a lu = ^m f e aflf3u Sj ia q p (A.7) 

and q = pi + P2 is the four-momentum of the virtual photon. Taking into account 
formula (|A.5|) , we get 



a. 



e atl p„V2n?q p . (A.8) 



^"3 X! 

Here P2 — £2-Pv^2^2 and 711(2) = 4=(1, 0, 0, ±1). Substituting the expression (|A.8| ) 



into the antisymmetric tensor ( |A.6[ ), and recalling the kinematic relations deduced in 



sect. 3, A* turns out to coincide with formula 
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Appendix B 

Here we derive the expression of the spin correlation matrix of two quarks, more- 
over we study in detail the coupling between a transversely polarized quark 
and a linearly polarized gluon. 

B.l - Spin correlation matrix 

The spin correlation matrix is defined as 

K T ' T '(p,p')=u T (p)u T ,(p'), (B.l) 

where Ux(p) is the Dirac spinor of a quark with four-momentum p and third spin 
component T. We consider a nonvanishing quark mass. First of all, we take both 
quarks at rest, secondly we make a Lorentz boost, thirdly we change the momentum 
of one of the two fermions; lastly we treat the limiting case of a negligible mass. 

B.l.l - Both quarks at rest 

Choosing, as usual, the z-axis as the quantization axis, we have 
2 

K o' T = Yl xjxj^ij = a 11 (l + o- 3 ) + a 12 (cr 1 + icr2) + a2i(o-i-icr2) + a22(l-cr 3 ), (B.2) 

where x T is the usual spinor. For pure spin states, such as those we are interseted in, 
some coefficients are nonvanishing: for T — T' a i2 and a 2 i vanish, while for T ^ T' 
On = a 2 2 = 0. In particular we have 

Kq ,t = m(l + 2T<t 3 ), k^~ t = m(a 1 + 2iTa 2 ), (B.3) 

having choosen for the Dirac spinor the normalization ut{p)ut{p) = 2m, where m is 
the rest mass of the quark. The previous expressions can be generalized by introducing 
three mutually orthogonal unit vectors, i, j and k. Choosing the quantization axis 
along k, we obtain 

a 3 -> Gik\ <7i + 2iTa 2 -> a t b T , b T = i + 2iTj. (B.4) 
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In terms of the Dirac matrices the spin correlation matrix at rest can be written as 



kY = ^ (1 ± 7 o)(l - 2T 75 7,n 4'" T = -y (1 ± 7o)75 (B.5) 

where m is the fermion rest mass and the ± sign in front of 70 refers to 
the (anti-)quark. 

B.1.2 - Quarks with equal momenta 

The previous formulae can be easily extended to the case of two quarks with equal 
momenta. We get, similarly to the spin density matrix, 



J p T = l -{i>±m){l + 2T 1 4), K T V >~ T = i#±m)7 5 0r, (B.(>) 



where, in the rest frame of the fermions, k = (0, k) and = (0, by). 

It is worth noticing that, if the two quarks have different third spin components 
and momenta not parallel to the quantization axis, the final result (second formula 
( |B.6| )) depends on the choice of the x and y-axis. We shall solve this ambiguity 



below, showing that the y — z plane must be fixed in such a way to contain the quark 
momentum. 

B.1.3 - Quarks with different collinear momenta 

Consider two quarks with different momenta, in a frame where they are collinear. 
In this case the correlation matrix reads 

K T > T '( P ,p') = U(p,p')K T /'. (B.7) 

Here p = (E, p) and p' = (E', p') are the quark four-momenta in the frame considered; 
moreover U(p,p') is the transformation matrix corresponding to the boost which 
changes p' to p, i. e., 

U(p, p') = -^== (l + j^j7o7**) , (B.8) 



where 

. Ip'I - IpI 



E' + E 
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(B.9) 



By using the Dirac equation, eq. ( |B.7| ) yields 



T,T'f /N 1 ft | + E'TjQm \ rprpl 

^' p ) = 7rrf l 1+t | P /| )V • (b.io) 



B.1.4- Limit of zero quark mass 

This limit is trivial in the case of X" = T: 

K^fojO = 2^^^ (1 + 2TT5 ^ (B ' 11} 
On the contrary, as we have observed above, if T 1 ^ T, and if the quark momenta 
p and p' are not parallel to the quantization axis, i. e. to k, we are faced with an 
ambiguity. In particular, let us consider the case when p is orthogonal to k, which is 
of interest to us: the correlation matrix is different according as to whether we take 
i or j along p. We show that the latter choice is correct. 
Take j along p. Formula ( |B.10| ) yields 

K T >- T (p,p') = ^=L= (l + t^p) (/±m) 75 (2zT/ c m- 1 +f), (B.12) 

where p' c = (|p'|, 0, 0, £") and I = (0, i). In the limit of m — > 0, taking into account 
that p"p" c = 0(m 2 /p' 2 ), we get 

K T '- T (p,p') = ^-^=/ 75 (T2zT + J/). (B.13) 

In the specific case that we consider in sect. 4, we identify the four-vector I with S, 
which is given by eq. (f|). In order to check our choice, consider two massless quarks 
with equal momenta. Taking the y-axis along the quark momentum, we have 

u 1 ' 2 = l=(u R + u L ), u- 1 ' 2 = ^={u R - u L ), (B.14) 

where ur and ul are respectively the spinors of the right-handed and of the left- 



handed quark. Substituting formulae (B.14Q into eq. (p.lj ), we get 

Kp~ T = --\2T(u R u R - u L u L ) + u L u R - u R u L }. (B.15) 
The usual formulae for the density matrix yield UrUr — Ujju,^ = ±^75. On the 



other hand, applying the second eq. ( p.6| ), we get ulUr — urUl = ifl^- Therefore 
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eq. ( |B.15| ) turns out to coincide with eq. ( |B.13| ) for t = 0. Now return to eq. ( |B.10| ) 



and take i (instead of j) along p. This choice yields a result which differs from eq. 
(|B.13|) by a factor 2iT, and it does not satisfy the self-consistency test that we have 



just illustrated. 

Notice that in the massless case t = {x' — x)/(x' + x), where x and x' are the 
light cone fractional momenta of the quarks; therefore the expression ( [B.13|) is covari- 
ant. We do not take into account the quark transverse momentum - which involves 
a Melosh-Wigner rotation -, however this is not requested in the twist-three approx- 
imation, which we are considering. 

Lastly it is straightforward to see that in eqs. (p.llj) and ( |B.13|) we can interchange 
p <-> p'. This property we exploit in subsect. 4.1. 

B.2 - Linearly polarized gluons 

Now we consider a gluon emitted by one of the two protons (say proton 2) and 
interacting perturbatively with the active quark of the other proton. In this connec- 
tion we observe that linearly polarized gluons are naturally coupled to transversely 
polarized quarks. In the reference frame defined above, for a gluon travelling along 



the y-axis, we have|49|, 50 



G% = -^(G R + G>t), G^ = -j=(G£-G£), (B.16) 

where G% and G^ z denote gluon states polarized respectively along the x- and along 
the z-axis and Gr and Gl are right-handed and left-handed gluons. Colour indices 
have been omitted for the sake of simplicity. The perturbative coupling of a massless 
quark with a gluon is of the type G^Hrj^ur or G^ulj^ul. On the other hand we 
have to do with quark states like ( |B.14j) , which couple with gluon states of the type 
( |B.16| ). Such couplings read 



(u Tl 7 ^ Tl )G£ = ^(G R u Rl ^ R + G>iu Ll ^u L ), (B.17) 

IT 

{u T ^u T ^)G» = -^(G R u Rl ,u R + G»u Ll ,u L ). (B.18) 
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Furthermore the gluon is related to the active quarks of protons 1 and 2 by angular 
momentum conservation, which implies that we have to take into account the factor 
F = sgn < 1/2, Ti; l,T g \l/2,T{ > = 27\, where T g = T' 2 - T 2 and T[ = T g + T\. 
Angular momentum conservation implies as well that for T' 2 = —T 2 we set T 2 = T\. 

Therefore the gluon polarization four-vector is e = 2T 2 / for T' 2 = T 2 and e = -k 
5ti,t 2 f° r T% = — ^2> where, in the reference frame defined in subsubsect. A. 1.4, I = 
(0,i) and k = (0,k). 



Appendix C 

Here we write the hadronic tensor according to quantum field theory, up to first 
order in g, adopting the axial gauge A a ■ n 2 = 0, with a = 1 .... 8: 

H G J = Hi ^V-^[r^(y)f^(x,y) + (l - 2)][r S(x-y)^r^', (C.l) 
where 



f%,(x,y) = <P0 p , i {U)D j i{y)^ j {x)\P l >, (C.2) 
D^y) = idJv+gAKy)^, (C.3) 

v l Uv) = <pM t mUy)\Pj>, (C4) 

x*-v) - l^^r^- ^ 

Two observations are in order. First of all in the above expression we have omitted 
the chronological products, according to the considerations by various authors |56| (see 
also[0). Secondly, the covariant derivative in the expression ( |C.2| ) implies gauge in- 
variance of the tensor ( |C.1|) . However we have omitted the gauge "link" operators [p2 



in the expressions (|C.2|) and (|C.4|) , since in the gauge adopted these operators can 
be made unity. All this simplifies a lot the calculations and immediately allows to 
establish symmetry properties of the correlation functions (see subsect. |4.2| ). 

We consider the usual Fourier expansions of the of the massless fermion and vector 
boson fields, i.e., 
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1 1/2 3. ^3 

M x ) = q^JIJZ J 2 i/ a|p |i/ a [4,i(p)«a(p)x«e qw 

+4!,(p)^(p)x,e-H, (C.6) 

^•a(^) = E/ ^KWCfkle'Hc,.]. (C.7) 

Here e m (m = 1, 2) are the two different polarization four- vectors of the physical 
gluons: e 1 = S (i. e., along the proton spin) and e 2 = S (i. e., orthogonal to 
the proton spin and to the proton momentum). a c m and are the creation and 
destruction operators of linearly polarized gluons. 

Substituting expressions ( |C.6| ) and ( |C.7| ) into eq. ( |C.1| ), we obtain a tensor of the 
type 

= H$£ + g(H^ a - H^ b ). (C.8) 

Here, recalling the definition 

£ < Pilc^JpiHjPi) \Pi >= qIt^upL), (C.9) 

8=1 

and an analogous one for antiquark densities, we find Hq^I = (see eq. (|22"D in 
the text). Moreover we can identify H^l' 11 ^ with H'^ a ^ (cfr. eq. (|40|)), provided we 
normalize the quark correlation functions according to eq. ( 1551) and assume a similar 
formula for the antiquark correlation functions. We stress that this normalization 
guarantees gauge invariance. 
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